Photo-assisted heat current and Peltier coefficient in a metal/dot /metal junction 
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The photo-assisted heat current through a metal/dot /metal junction and its associated Peltier 
coefficient are computed in the framework of the time-dependent out-of-equilibrium Keldysh for- 
malism in the presence of a dot energy modulation. When the frequency of the modulation is much 
larger than the amplitude of the modulation, the heat current follows the sinusoidal time evolution 
of the dot energy. This is no longer the case when the modulation frequency becomes of the order 
or smaller than the amplitude of the modulation. To characterize this non sinusoidal behavior, we 
have calculated the harmonics of the photo-assisted heat current. The zero-order harmonic can 
be expressed as an infinite sum of dc heat currents associated to a dot with shifted energies. It 
exhibits a devil staircase profile with non horizontal steps whereas it is established that the steps 
are horizontal for the zero-order harmonic of the photo-assisted electric current. This particularity 
is related to the fact that the dot heat is not a conserved quantity due to energy dissipation within 
the tunnel barriers. 
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I. INTRODUCTION 

The photo-assisted electric current through an insu- 
lating barrier is the topic of a large number of stud- 
ies both theoretically[l]-|4| and experimentally in a 
wide range of systems, such as for example: quantum 
dots, normal or superconducting tunnel junctions, and 
Luttinger liquids. In the presence of a voltage modula- 
tion, the electrons can emit or absorb photons when they 
travel across the barrier and the resulting current is thus 
a superposition of dc currents. When the current-voltage 
characteristic is non-linear, it leads to a specific behavior 
for the current. @ 

This is only recently that the study of photo-assisted 
heat current has emerged. On the theoretical side, the 
heat flow generated by an adiabatic quantum p ump 
through a mesoscopic sample has been studied @, [l0| . 
as well as the photo-assisted heat flow in a nor- 
mal metal/superconductor i unction [ill]. Non-adiabatic 
pumping heat in an asymmetric double quantum dot has 
also been considered [l 2] . Even more recently, the chiral 
heat transport in driven quantum Hall edges states fl3j 
and the microwave-mediated heat transport and ther- 
moelectric effect in a q uantum dot in the presence of 
Coulomb interaction [lj] have been calculated. On the 
experimental side, measurements of Seebeck voltage has 
been performed in magnetic tunnel junctions in the pres- 
ence of a frequency modulated laser used to heat up the 
device [lEj. 

The formalism that has been developed [l6| to calcu- 
late the time-dependent heat current with the help of 
out-of-equilibrium Green's functions allows to treat any 
type of time-dependent voltages, in particular a modu- 
lated one. It is therefore the appropriate formalism that 
has to be used when one considers time-dependent exci- 
tation in the quantum regime. It enables not only to cal- 



culate the photo-assisted electric and heat currents but 
also, as done here for the first time, the photo-assisted 
Peltier coefficient since this last quantity is defined as the 
ratio between the heat current and the electric current 
in the absence of temperature gradient. However, this 
formalism does not make possible the calculation of the 
photo-assisted Seebeck coefficient. Indeed, whereas this 
quantity is generally measured in an open circuit, it is 
not straightforward in such a calculation to ensure the 
cancellation of the electric current, except in the linear 
regime [171] which is not the one we consider here. 

Even if several theoretical works are devoted to the 
photo- assisted heat current in quantum dotsfTH, [IJ, Il8r - 
[20 ] , they are restricted to the calculation of its time aver- 
age value which corresponds to the zero-order harmonic. 
A detailed study of higher order harmonics is still miss- 
ing. The main objective of the present work is to fill 
this lack. This will allow us to know if the heat cur- 
rent through the junction will follow the time evolution 
of the modulated gate voltage, i.e., if it is sinusoidal, or 
if it has a more complicated profile. A careful study of 
the heat current and Peltier coefficient according to the 
amplitude and the frequency of the modulation, as done 
here, is thus needed. 

The paper is organized as follows: in Sec. II, we expose 
the model we have used to calculate the time-dependent 
heat current. In Sees. Ill and IV, we discuss the photo- 
assisted heat current and the photo-assisted Peltier co- 
efficient. Next, in Sec. IV, we study in more details the 
harmonics of the ac heat current. We conclude briefly in 
Sec. V. 



II. MODEL 

We consider a single level non-interacting quantum dot 
connected to left (L) and right (R) reservoirs with chem- 
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ical potentials (J.l,r at temperatures Tl,,r (see the upper 
panel of Fig. [T]). We define as well the source-drain volt- 
age: eV = fiL — M-Ri ancl t ne average temperature of the 
reservoirs: T = (Tl + T R )/2. The energy level of the 
quantum dot can be modulated in time by tuning the 
gate voltage, we denote it £dot(*)- We use the following 
Hamiltonian to describe this system: 
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k£L,R 



£dot 
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(1) 



where £keL,R. is the energy band of the reservoir L or R. 
The notations c\. (eft) and Ck (d) refer to the creation and 
annihilation operators associated to the reservoirs (dot). 
VkeL,R is the hopping amplitude between the reservoir L 
or R and the dot. 
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FIG. 1: (Upper panel) Schematic picture of the left (L) and 
right (R) reservoirs connected to a single level quantum dot 
with a modulated gate voltage. Notice that the definition we 
take for the currents means that their directions are oriented 
from each of the reservoirs to the dot. (Lower panel) Picture 
of the equivalent energy levels of the dot. The heat current is 
sinusoidal in the case where only the energy levels represented 
by a thick line contribute (i.e., the energy levels eo and £o±w). 



The time-dependent heat current through the left (L) 
or right (R) reservoirs, l\ R = —Ql,r, is related to the 
energy current, if R = —El,r, and to the electric cur- 



rent, I? p = —cNl,r, by the identity: 



VL,R 



which is obtained from the thermodynamic relation: 
dQL,R = dE L ,R - Li LtR dN LiR , 



(2) 



(3) 



since the reservoirs are at equilibrium. The quantities 
Ql,r, El R and Nl r correspond respectively to the 
heat, the energy and the number of particles in the 
left (L) or right (R) reservoirs. Thus, by calculating the 



electric and energy currents, one can deduce the heat 
current with the help of Eq. @. To determine these 
quantities, we use the fact that ihNL,R — [Nl,r, H] and 
iHE LtR = [E LiR ,H] with: 
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Notice that the definitions of the energies of the reser- 
voirs, El iR , given above implicitly assume that we con- 
sider only the electronic contribution to the energy cur- 
rent (i.e., we neglect possible phononic contribution). 

Within this model and with the help of Kelysh out- 
of-equilibrium Green's functions, the expression of the 
time-dependent heat current through the reservoir p in 
the wide band limit reads as: 161 
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(e - Hp)fp(e)lm{A(s, t)}ds 
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(e - ^ P )f P '(e)\A(e,t)\ 2 de 



(6) 



where f p is the Fermi-Dirac distribution function of the 
reservoir p. The enlargement of the dot energy level 
due to its coupling with the reservoir p is defined as 
T p = 2n p p \Vp\ 2 , where p p is the density of states of the 
reservoir p, and V p = Vfcg p . In the wide band limit, it 
is assumed to be energy independent. Notice that this 
is the energy current contribution, I*j — —E p , which is 
responsible of the presence of the e term in the heat cur- 
rent. 

The quantity A(e, t) in Eq. © is the spectral function 
defined as:[2l| 



A(e,t) 



G^tMV^'^^dh 



(7) 



The retarded Green's function G r is the one of the dot 
connected to the reservoirs which is given here by: 

G r {t,t l )=9 r {tM)e T{tl - t)/2h , (8) 

where g T {t,t x ) = -*6(* - ti)e" l/t 'i dt2£ *°^l h i s the re- 
tarded Green's function of the isolated dot, r = Tl +T r , 
and <d is the Heaviside function. 

In the present work, we consider the following time 
modulation for the dot energy level: 

£dot(*) = £o + £i cos(wt) . (9) 
In that case, the spectral function of Eq. (0) becomes: 



A{e,t) = 



E 



Jn( — 



(- 



E n {e)+i^ 



(10) 

where E n (e) = e — £o — nu, and J n is the Bessel's 
function of order n which appears in the expression of 
the spectral function through the identity e la;sm M) = 
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III. PHOTO-ASSISTED HEAT CURRENT 

In this section, we look at the time evolution of the 
heat current Ip(t) which has been obtained by insert- 
ing Eq. (|10[) in Eq. (j6|) and integrating over energy nu- 
merically We discuss its principal characteristics and 
compare them to the ones of the time-dependent electric 
current Ip(t) that has a similar form than Eq. (|6]) except 
the factors (e — fi p ) which are not present. |2l| 




FIG. 2: (Upper panel) Heat current and (lower panel) electric 
current as a function of uot. The parameters are: eo = 1, £i = 
0.5, Yl,r = 0.5, k B T L , T = 0.01 and eV = 1. The modulation 
frequency is: hw = 10 (solid red lines) and hw = 0.5 (dotted 
green lines). The black straight lines correspond to the dc 
currents in the absence of modulation (i.e., ei = 0). The 
thick (thin) lines correspond to the photo-assisted currents in 
the left L (right R) reservoir. The unit for energies is F. 



It is important to notice that there are several char- 
acteristic energies in this system that are related to: the 
temperatures, /cbIl^, the source-drain voltage, eV, the 
modulation frequency, Huj, the amplitude of the mod- 
ulation, £]_, the dot energy level, So, and the coupling 
strength between the dot and the reservoirs, T. The heat 
current depends on the relative values of all these ener- 
gies. Here, we focus on the change in its behavior accord- 
ing to the value of the ratio ei/fku that appears in the 
argument of the Bessel's functions present in the expres- 
sion of the heat current. Fig. [2] shows the time evolution 
of the heat and electric currents, at fixed E\ and temper- 
atures HbTl r much smaller than all the other energies, 
for several values of the modulation frequency u during 
a full period. We see that when &j > £[, the left and 
right currents are both sinusoidal (see the solid red lines) 
and oscillate almost in phase around their correspond- 
ing dc values (see the black straight lines) whereas when 



hu> < Ex, the signals are much more complex (see the dot- 
ted green lines). Indeed, when the ratio Si/Huj is smaller 
than one, the dominant contributions in the heat current 
come from the terms associated with n, m G { — 1, 0, 1} in 
the sums that appear in the spectral function of Eq. (|T0|) 
because the value of the product of Bessel's functions, 
Jn(x)J m (x), with argument x smaller than one decreases 
quickly with increasing n and m (see Fig. 13]). As a con- 
sequence, we have: 



Th(0) 

V 



2Re{/. 



h(l) Aojt 



} 



(11) 



where Ip^ is the zero-order harmonic (or average value) 
of the time-dependent heat current, and Ip' its first 
order harmonic (see Sec. V). This result explains the si- 
nusoidal behavior we have obtained in the limit tuo 3> £\ 
(i.e., the red solid lines of Fig. [5]). On the contrary, when 
the ratio Si/fkJ is close or higher than one, a large num- 
ber of terms in the sums over n and m in Eq. (110[) con- 
tribute and Ip (t) is a superposition of several harmonics 
(see Sec. V) which lead to a non sinusoidal signal (i.e., 
the green dotted lines of Fig. [2]). 




FIG. 3: Absolute value of the product J n (x)J m (x) that ap- 
pears in the spectral function of Eq. (|10l) as a function of 
argument x for different values of n and m. When x < 1, 
this product decreases quickly when n or m increases. On the 
contrary, when x > 1, there is no dominant terms in the sums 
over n and m. 



IV. PHOTO-ASSISTED PELTIER 
COEFFICIENT 

We have all the ingredients to calculate the time evo- 
lution of the Peltier coefficient of the junction which is 
defined such as: [la 



n(t) 



(12) 



By looking on Fig. [2l we notice that whereas the left 
and right heat currents can be equal at some specific 
times (compare the thick and thin lines in the upper 
panel of Fig. [5]), this is not the case for the left and right 
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electric currents which take distinct values at any time 
(compare the thick and thin lines in the lower panel of 
Fig. [2]). This property is general and not limited to the 
parameters we have chosen for plotting Fig. [5J Thus, 
we can conclude that the Peltier coefficient of the junc- 
tion will never be divergent. However, it vanish every 
times that the left and right heat currents are equal, i.e. 

Figure @] shows the time evolution of the Peltier co- 
efficient at low temperature (i.e., at feeTr,^ <C V). As 
it was the case for the heat and electric currents, the 
Peltier coefficient is sinusoidal in the limit huj S\ (see 
the solid red line): it follows the time evolution of the 
imposed gate voltage modulation. When the modulation 
frequency is reduced and becomes of the order of ampli- 
tude of the modulation fkj ~ Si, the Peltier coefficient is 
no longer sinusoidal (see the dotted green line). 
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FIG. 4: Peltier coefficient as a function of uit. The parameters 
are the same to those of Fig. [2] eo = 1, ei = 0.5, Tl,r = 0.5, 
fcsTi^ = 0.01 and eV = 1. The modulation frequency is: 
liw — 10 (solid red line), and hw = 0.5 (dotted green line). 
The black straight line corresponds to the Peltier coefficient in 
the absence of modulation (i.e., ei = 0). The unit for energies 
is r. 



Another quantity which characterizes the thermoelec- 
tricity is the Seebeck coefficient defined as the ratio be- 
tween the voltage gradient and the temperature gradient 
in an open circuit: S = AV/ AT\j^ = q. However, this 
definition applies only in the linear response regime, as 
well as the direct relation between Peltier and Seebeck 
coefficients: n = ST. Since here we are not necessar- 
ily in the linear response regime, the calculation of the 
photo-assisted Seebeck coefficient is a more difficult task 
which goes beyond the scope of the present study. 



HARMONICS OF THE HEAT CURRENT 



We turn now our attention to the harmonics, I 



h(N) 
L,R 



of 



the left and right heat currents which are defined through 
the relation: 



oo 
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AT=— oo 



(13) 



The fact that l£ R (t) is a real quantity imposes that 



I 



h(N) 



, /; and 1^ r are complex conjugates. 
Reporting Eq. (fTQ"j) in Eqs. © and (TlU)) . and perform- 
ing the integration over time, we obtain: 



T h(N) 
l L,R 



L L,R. 

~2h~ 



± n— — oo 



f- 



rfi,i(/i,i?(e) - T 2if LM {e)E n±N (e 

£ — V>L,R 



(E n (e)^ir/2)(E n±N (e)±iT/2) 



>de 



(14) 



where we have used the identity: ~^2 n J n + rn {x)J n (y) = 

Jm(x ~ y). 

In the following, we study separately the zero-order 
harmonic of the heat current and the its harmonics of 
higher orders. 



A. Harmonic of order zero 

The zero-order (N — 0) harmonic of the heat cur- 
rent corresponds to the time average heat current. From 
Eq. (fl~4)) . we obtain: 



T h(0) 
L,R 



n—~oo 



^(^VSW^), (15) 



\huj ) 



where is the dc heat current defined as: 



iTr\z) 



(g - ML,fl)(/L,fl(£) - fR,U £ )) 

(e-2) 2 +r 2 /4 



ds 



(16) 



The description of Eq. flT5l) is the following: the zero- 
order harmonic of the ac heat current is equal to an infi- 
nite sum over n of dc heat currents associated to shifted 
dot energies: Eq + nuj, times the square of the Bessel's 
function of order n (see the lower panel of Fig. Q] for a 
schematic representation of those shifted energy levels). 
A similar kind of relation linking the zero-order harmonic 



of the electric current, I, 



.(0) 



and dc electric currents, 



Ip (dc) {e + nuj), holds 



7- e (°) _ 

l L,R ~ 



where: 



Ti — — OO 



4 (£)«»(« 



UMe) - Ir,l{ £ ) 
(e-z) 2 +T 2 /4 



ds 



(17) 



(18) 



Even if the expressions of the zero-order harmonics 
of the heat and electric currents look similar, i.e. they 
contain an infinite sum of dc contributions, they differ 
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on an essential point: whereas I L 



e(0) 



r e(0) 



we have 



I 



Ho) 



i ^= ~Ir because of the (e — ^l.r) factor which 



is present in Eq. (TTBl) but not in Eq. (|TB]). This leads 
to rather different profile for the zero-order harmonics of 
the heat and electric currents as detailed below. 
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FIG. 5: Zero-order harmonic of the left heat current J^' ' 
(upper panel) and right heat current /]j (lower panel) as a 
function of eo. The temperatures are /cbTl,t = 0.01 (solid 
black lines), ksT^.T = 1 (dashed blue lines) and fcsTL.r = 2 
(dotted orange lines). The other parameters are ei = 10, 
eV — —10, and frw — 10. The unit for energies is Y. 



Figure [5] shows the zero-order harmonic of the left and 
right heat currents at fixed frequency modulation, hui = 
10r and voltage, eV = — 10r, as a function of £o- We 
observe a jump each time that the energy (e ± V/2) 
is an integer multiple of Huj, here: £o = 5T, Eq = 15T, 
etc... Between these points, the left (right) heat current 
slowly decreases (increases) contrary to what is observed 
for electric current which takes a constant value between 
the jumps (see Fig. [S]). The staircase behavior of the 
zero-order harmonic of the electric current is a standard 
result 3 in the regime corresponding to the one of Figs. [5] 
and [5] (i.e., T <C {Hoj, e\, eV}). It comes from the fact 
that the dot occupation is a conserved quantity in the 
stationary regime, contrary to the dot heat which is not 
conserved, even in the stationary regime, because of the 
energy dissipation within the tunnel barriers (l6i Il7j . 

The width of the jumps depends both on the temper- 
atures and on the coupling strength between the quan- 
tum dot and the reservoirs (i.e., on the value of T). The 
higher the temperatures, the larger the jump width (com- 
pare the solid black, dashed blue and dotted orange lines 
in Figs. [5] and H]). The increasing of the coupling T 
enhances also the larger of the jump width and when 
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Zero-order harmonic of the electric current T 
as a function of eo. The parameters and the legends 



are the same to those of Fig. [5] The unit for energies is V. 



r ~ {fvjj, £i, eV}, the staircase structure is not more vis- 
ible (not shown). 

B. Harmonic of order N 

With the help of Eq. (HU, we have calculated numer- 
ically the harmonics of the heat currents. In Fig. [7] is 
plotted the module of the harmonics of order 0, 1 and 2 
as a function of the modulation frequency. The module 
of the higher order harmonics of the heat currents, i.e. 



T h(N>2) 



, are not shown on the graphs because they are 

fe(2). 



much smaller in amplitude in comparison to \I L R 
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FIG. 7: Amplitude of the zero-order (solid black lines), first 
order (dashed red lines) and second order (dotted purple lines) 
harmonics of the left heat current (upper panel) and right heat 
current (lower panel). The parameters are ksTL,R = 0.01, 
eo = 1, eV — 1, and ei = 0.5. The unit for energies is V. 

A crucial information that we extract from Fig. [7] is 
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the fact that when Hu> < e%, the modules of first and 
second harmonics are of the same order in magnitude: 

l-^iitl ~ l-^ifll' whereas when hw 3> £1, the first har- 
monic dominates over the second harmonic (compare the 
dashed red lines to the dotted purple lines). These last 
result justifies the approximation that we have made in 
order to write Eq. and explain the sinusoidal behav- 
ior of the time-dependent heat current obtained in that 
regime. 

VI. CONCLUSION 

We have highlighted several interesting features in the 
photo-assisted heat current through a metal/dot /metal 



junction: (i) its zero-order harmonic exhibits a devil 
staircase with non horizontal steps which result from the 
fact that the dot heat in not a conserved quantity even 
in the stationary regime because of energy dissipation 
within tunnel barriers; (ii) the time evolution of the heat 
current is non sinusoidal when the modulation frequency 
of the gate voltage is of the order of the amplitude of 
the modulation; and (iii) the time evolution of the as- 
sociated Peltier coefficient is strongly dependent of the 
energy scales which characterize the junction. 
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